arXiv:1506.01988vl [physics.flu-dyn] 18 Dec 2014 


Papers in Physics, vol. 6, art. 060006 (2014) 


Received: 20 March 2014, Accepted: 7 August 2014 
Edited by: A. Marti 

Licence: Creative Commons Attribution 3.0 
DOI: http://dx.doi.org/10.4279/PIP.060006 


www. papersinphy sics.org 



ISSN 1852-4249 


Influence of surface tension on two fluids shearing instability 
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Using extended Layzer’s potential flow model, we investigate the effects of surface tension 
on the growth of the bubble and spike in combined Rayleigh-Taylor and Kelvin-Helmholtz 
instability. The nonlinear asymptotic solutions are obtained analytically for the velocity 
and curvature of the bubble and spike tip. We find that the surface tension decreases 
the velocity but does not affect the curvature, provided surface tension is greater than a 
critical value. For a certain condition, we observe that surface tension stabilizes the motion. 
Any perturbation, whatever its magnitude, results stable with nonlinear oscillations. The 
nonlinear oscillations depend on surface tension and relative velocity shear of the two 
fluids. 


I. Introduction 

When two different density fluids are divided by 
an interface, the interface becomes unstable with 
exponential growth under the action of a constant 
acceleration acting in the direction perpendicular 
to the interface from the heavier to lighter fluid or 
under the action of relative velocity shear of two 
fluids. These two types of instabilities are known 
as Rayleigh-Taylor and Kelvin-Helmholtz instabil¬ 
ities, respectively. Temporal development of the 
nonlinear structure of the interface consequent to 
Rayleigh-Taylor or Kelvin-Helmholtz instability is 
currently a topic of interest both from theoretical 
and experimental points of view. The nonlinear 
structure is called a bubble if the lighter fluid pen¬ 
etrates across the unperturbed interface into the 
heavier fluid and it is called a spike if the opposite 
takes place. The instabilities arise in connection 
with a wide range of problems ranging from direct 
or indirect laser driven experiments in the abla- 
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tion region at compression front during the pro¬ 
cess of inertial confinement fusion mia to mix¬ 
ing of plasmas in space plasma systems, such as 
boundary of planetary magnetosphere, solar wind 
and cluster of galaxies [3]. In high energy den¬ 
sity physics (HEDP), formation of supernova rem¬ 
nant or formation of astrophysical jets BHEI are 
also seen in these types of instabilities. In high 
energy density plasma experiments using Omega 
laser [9], Kelvin-Helmholtz instability growth has 
recently been observed . 

There are several methods to describe the non¬ 
linear structure of the interface of two constant 
density fluids under potential theory and the as¬ 
sociated nonlinear dynamics has been studied by 

many authors HDHI31- Layzer m described the 

formation of the structure using an expansion near 
the tip of the bubble or the spike up to second 
order in the transverse coordinates in two dimen¬ 
sional motion and this approach was extended in 
Ref. [H] for Kelvin-Helmholtz instability. It is 
well known m that the surface tension reduces 
the linear Rayleigh-Taylor Growth rate. The low¬ 
ering in the growth rate is seen to increase with 
increase in the wave number k up to a critical wave 
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number k c = \J , where T denotes surface 

tension, and ph and pi are the densities of the heav¬ 
ier and lighter fluids, respectively. The same effect 
has been described by Mikaelian m for Rayleigh 
Taylor instability in finite thickness and Sung-Ik 
Sohn [17] described the effect using the Layzer non¬ 
linear potential model. The nonlinear theory influ¬ 
ence of surface tension was elaborately studied by 
Pullin pTS] and Gamier et al. IS] using numerical 
methods. 

The present paper addresses to the problem of 
the time development of the nonlinear interfacial 
structure caused by combined Rayleigh-Taylor and 
Kelvin-Helmholtz instability in presence of surface 
tension. It is shown that the growth rate of the 
instabilities is affected by the surface tension. The 
growth rate of the tip of the bubble or spike are sig¬ 
nificantly reduced due to the surface tension. We 
observed an oscillatory stabilization of the inter¬ 
face for large surface tension. This oscillation de¬ 
pends on the relative velocity shear also. Section II 
deals with the basic hydrodynamical equations to¬ 
gether with the geometry involved. Here we assume 
that the fluids are inviscid and the motion is irrota- 
tional. The investigation of the nonlinear aspect of 
the structure of the two fluids interface is facilitated 
by Bernoull’s equation together with the pressure 
balance equation at the interface. The long time 
asymptotic behavior of the bubble and spike tip for 
combined Rayleigh Taylor and Kelvin Helmholtz 
instabilities is derived in section III.A and III.B, 
respectively. We have also discussed the character¬ 
istics of the tip of the bubble and the spike derived 
analytically and numerically. Finally, we have con¬ 
cluded the results in section IV. 

II. Basic mathematical model 

We have considered two incompressible fluids sep¬ 
arated by an interface located at y = 0 in a two- 
dimensional x — y plane, where x axis lying nor¬ 
mal to the unperturbed fluid interface. The fluid 
with density ph is assumed to overlie the fluid with 
density pi and gravity is taken along negative y- 
axis. In the following discussion, we shall denote 
the properties of the fluid above the interface by 
the subscript h and below the interface by the sub¬ 
script l. After perturbation, the nonlinear interface 
is assumed to take up a parabolic shape, given by 


y = r)(x,t) = 770 (f) + p 2 {t){x - 771 (f )) 2 ( 1 ) 

The perturbed interface forms a bubble or spike 
according to 770 (f) > 0 , 772 (f) < 0 or 770 (f) < 0 , 
b 2(0 > 0. Functions 770 (f) and 771 (f) are related 
to the position of the tip of the bubble from the 
unperturbed interface, i.e, at time t the position of 
the bubble tip is ( 771 (f), 770 (f)) and 172 (f) is related 
to the bubble curvature. 

In our previous works |14 [i20H^5] . we have con¬ 
sidered 771 (f) = 0 due to the absence of velocity 
shear parallel to the unperturbed interface. How¬ 
ever, in presence of streaming motion of the fluids, 
the tip of the bubble moves parallel to unperturbed 
interface with velocity 771 (f). 

According to the extended Layzer model mm 
it™ , the velocity potentials describing the motion 
for the upper (heavier) and lower (lighter) fluids are 
assumed to be given by 

4>h{x, y, f) = ai(f) cos {k(x - p 1 (t))e~ k(v ~ Vo(t)) 

+ 02 (f) sin (k(x — pi{t))e~ k< ' v ~ r,o( ' t ' ,Sl 

- xU h (2) 

c/)i(x,y,t) = b 0 (t)y 

+ 61 (f) cos ( k(x — 77 i(f))e fe ^ _?)0 ^^ 

+ 62 (f) sin ( k(x — r]i{t))e k ^ v ~ r,0 ^ S> 

- xUi ( 3 ) 

where Uh and Ui are streaming velocities of up¬ 

per and lower fluids, respectively, and k is the per¬ 
turbed wave number. 

The evolution of the interface y = rj(x,t) can 
be determined by the kinematical and dynamical 
boundary conditions. The kinematical boundary 
conditions are 

dp _ djd^h _ /.n 

dt dx dx dy 

._d(h_d<k ,,, 

dx dx dx dy dy 

and the dynamical boundary condition (first inte¬ 
gral of the momentum equation) is of the form 
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— Ph(l )—- k ~Ph(l){v <Ph(l)) + Ph(l)9U 

= ~Ph(l) + fh(l){t) (6) 

The pressure boundary condition at two fluid in¬ 
terface including surface tension H3H23 is 

T 

Ph ~ Pl= R ( 7 ) 

where T is the surface tension and R is the radius 
of curvature. 

Plugging the condition (7) at the interface y = 
r)(x, t ) in Eq. ( 6 ), we obtain the following equation. 

T 

+g(ph- pi)y = -R + fh~ fi (8) 


k 2 b\ _ 66 + 1 

Vkg~ 66-1 


( 14 ) 


fe 2 b 2 _ (26 + 1)6 ~ 2 UVh ~ Vi) 

\fkp 2^-1 


dr 


^ 1 ( 6,6 a 

Di(&,r) (66-1) 

, 2(1 — 66(66 ~ 1) ( 
Oi( 6,6 \ 

^2(6,0 (66-i)4 2 

Or (6 >6 26(26-1) 
, 2(46-1)(66-1) 


1 - 126 


fc? 


A(6,r)(26-l) a 

x [(14 - 6) 2 6 - (14 - V))(26 +1)6] (16) 


and 


We have restricted our study near the peak of the 
perturbed structure where |fc(x— 77 i(t))| <C 1. Thus, 
we can neglect the terms of 0(\x — i]i\ l ) (i > 3) [13] . 
With this point of view, we have 

i = 2? ?2 (1 + 4 ? 7 |( a ;- ? 7 i ) 2 )' 5 

~ 2 ?? 2 (l — 6 ? 7 |(a; — rji) 2 ) (9) 

We substitute all the parameters rj 1 (f>h and <pi in the 
kinematic and dynamic boundary conditions repre¬ 
sented by Eqs. (4), (5), ( 8 ) and (9), and equate co¬ 
efficients of (x — ? 7 i)*,(i = 0 , 1 , 2 ) and neglect terms 
0{\x — rji\ l ) (i > 3). This yields the following equa¬ 
tions. 


d 6 

dr 


= 6 


( 10 ) 


^6 _ (26 — i ) r 66 

dr 26-D2(6>6 

6(66 + 1 ) 

20 2 (6,r)(66-l)(26-l) 

X [4(6, - 6)(46 - 1) - |(286 2 


(17) 

46 - 1 )] 


where r = 6 = krj 0 ; 6 = kr] i; 6 = ^; 

u = & = % r = t-M-, kl = fcato 

and V h( i) = kp are corresponding dimensionless 
quantities. The function -/Vi, 2 ( 6 i r ) and 01 , 2 ( 6 ) r ) 
are given by 


Ni( 6, r) = 36(1 - r)6 2 + 12(4 + r)6 + (7 - r); 

Oi( 6)6 = 12(r - 1)6 + 4(r - 1)6 -6 + 1) 

(18) 


<%2 Tr 6(26 + 1 ) 
d^ = Vt 2(T~ 


f=4< 6&+i >«‘ 

kbp _ 1266 

v^ _ 66-1 


and 

( 11 ) 

JV 2 (6, r) = 16(1 - r)6 3 + 12(1 + r)6 - (1 + r); 
(12) 0 2 (6,r) = 2(l-r)6 + (r + l) (19) 

The temporal development of the combined effect 
of Rayleigh-Taylor and Kelvin-Helmholtz instabil- 
ity is given by Eqs. (10)—(12), (16) and (17). 
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Figure 1: Bubble- variation of £i, £ 2 , £ 3 , £4 and 
£5 with r. Initial value £1 =0.1, £2 =0, £3 =-0.05, 
£4 = 0, and £5 = 0 with p h = 3, pi = 2, 14=0.5, 
Vi = 0.1, fi 2=0 (line), 0.5 (dot), 1 (dash), 3.9 (dash- 
dot). 


III. Numerical results and discus¬ 
sions 

i. Effect of surface tension on bubble 
growth 

In this section, we present the effect of surface 
tension on the nonlinear growth rate of the bub¬ 
ble tip for combined Rayleigh Taylor and Kelvin 
Helmholtz instability. To describe the dynamics 
of the bubble tip, it is essential to integrate Eqs. 
(10)-(12), (15) and (16) by numerical simulation. 
To obtain the initial conditions of the numerical 
integration, we assume that the initial interface 
is given by y = po(t = 0 )cos(kx). The expan¬ 
sion of the cosine function gives (£ 2 )initial = 0 and 

(£3 ^initial — y (Cl )initial ? where (£l )initial 1® the 

arbitrary initial amplitude. Since the perturbation 
starts from rest, we may often choose (£ 4 )initial = 
(£ 5 )initial = 0. The non-dimensionalized time de- 


Figure 2: Bubble- variation of £ 1 , £ 2 , £ 3 , £4 and 
£5 with r. Initial value £1 =0.1, £2 =0, £3 =-0.05, 
£4 = 0, and £ 5=0 with p h = 3, pi = 2, 14=0.5, 
Vi = 0.1, fr=10 (line), 15 (dot), 20 (dash). 


velopment plots of £ 1 , £ 2 , £ 3,£4 and £5 are shown in 
Figs. 1, 2 and 3. 

Before we describe the nature of the bubble 
tip, consider the asymptotic behavior of the tip. 
As r —> 00 , the asymptotic values of £ 3 , £4 and 
£5 for bubble are obtained by setting -p = 0 , 
= 0 and = 0. Note that, if k 2 < 
3 (l + i^(AP) 2 ) fc c 2 , where AV = 14 - V h 
the asymptotic values are 


[(6 


3 ) asymp\bubble 


1 

6 


( 20 ) 


K?4) 


asymp\bubble 


2 A 


3(1 + A) 



and 


( 21 ) 


5 l — A 
16 1 + A 


(AV) 2 


[(6 


5 JasympJ bubble 


= 0 


( 22 ) 
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Figure 3: Bubble- variation of £i, £ 2 , £ 3 , £4 and 
£5 with t. Initial value £1 =0.1, £2 =0, £3 =-0.05, 
£ 4 = 0, and £5 = 0 with p h = 3, pi = 2, p=20, 
V h =0, Vi = 0 (line), 14=0.5, V) = O.l(dot), 14=1, 
Vi = 0.1 (dash), 14=1.5, Vi = 0.1 (dash-dot). 

where, A = ph ~ p ‘ is the Atwood number. 

It is clear form Fig. 1 that surface tension sup¬ 
presses the velocity and growth of the bubble tip 
significantly, provided surface tension is larger than 
a critical threshold, T < T bubble , where 

T b c ubble = 3 (V - Pl ) + ^U(AP) 2 ) ± (23) 

Here the critical value depends on the magnitude of 
relative velocity shear of two fluids and the growth 
and velocity of the tip reduced if T < T bubble . 

When there is no tangential velocity difference 
(i.e., Vh = Vi) between the two fluids initially, the 
fluids are purely prone to the Rayleigh-Taylor in¬ 
stability and the critical value becomes 3 ^ Ph k2 Pl g • 
These results agree with the argument in Ref. CZI- 
In absence of surface tension, the asymptotic values 
coincide with the results obtained in our previous 
work El- 


Further, if T > T bubble , oscillatory state emerges 
even for r > 1. Figures 2 and 3 describe the os¬ 
cillatory state of the motion. The amplitude and 
the period of oscillation decrease monotonically for 
large surface tension (Fig. 2), while the ampli¬ 
tude of oscillation increases for large relative ve¬ 
locity shear (Fig. 3). In this respect, Figs. 2 and 
3 show that there always exists a self generated os¬ 
cillatory transverse velocity component (—£ 5 ) due 
to perturbation and this depends upon surface ten¬ 
sion as well as the relative velocity shear AV at 
the two fluids interface. For negative velocity shear 
(i.e, AV < 0), the self generated oscillatory trans¬ 
verse velocity of the bubble peak acts opposite to 
the direction of 14 and the amplitude of oscillation 
increases for large surface tension. 

If T = T bubble , equilibrium is attained, i.e, 

£3 = = £5 = 0 

when £3 = ~ and £4 = £5 = 0 (24) 

6 

and the equilibrium becomes unstable. This feature 
is shown with a dot-dash line in Fig. 1. Thus, the 
combined Rayleigh-Taylor and Kelvin Helmholtz 
instability is stabilized when 

t2>3 f 1 + i- a T (A ' / > 2 H 

v 16 Ph - Pi ) 

i.e., T > T bubble (25) 

while the instability however persists but with re¬ 
duced growth rate for 

P<3(l + l|-5!-(AC)be 

v 16 Ph - Pi ) 

i.e., T < T bubble (26) 

According to the condition (25), for ph = 3, 

pi = 2, 14 = 0.5 and V) = 0.1, the motion is stabi- 

.2 

lized when 4- > 3.9. These results are exhibited in 

Fig. 2, where > 3.9. For = 3.9, the growth 
rate of the instability is asymptotically diminished 
and becomes 0 (dash-dot line of Fig. 1). However, 
Fig. 1 shows the suppression of growth rate of the 

7,2 

instability due to surface tension, when p- < 3.9. 
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Figure 4: Spike- variation of £i, £ 2 , £ 3 , £4 and £5 
with r. Initial value £1 =-0.1, £2 =0, £3 =0.05, £4 = 
O^and £5 = 0 with p h = 3, pi = 2, 14=0.5, Vj = 0.1, 
p-=0 (line), 0.5 (dot), 1 (dash), 4.35 (dash-dot). 


Figure 5: Spike- variation of £ 1 , £2, £3, £4 and £5 
with r. Initial value £1 =-0.1, £2 =0, £3 =0.05, 
£4 = 0, and £ 5=0 with p h =3 , pi = 2, 14=0.5, 
Vi = 0.1, p-=10 (line), 15 (dot), 20 (dash). 


ii. Effect of surface tension on spike growth 

The temporal evolution of spike state is exhibited 
in Figs. 4, 5 and 6 ; the results follow from the 
numerical integration of Eqs. (10)—(12), (15) and 
(16) using the transformation £1 -A —£1, £3 —»• —£3, 
g —g, r —> £ and 14 V). The saturation 

curvature and velocity of the spike tip are given by 


[(^3') asymp] spike — g 

(27) 

[(£4 asymp] spike 

(28) 


Figure 4 describes that large surface tension sup¬ 
presses the growth rate of the spike tip, as well as 
the bubble. The nonlinear oscillation of the spike 
tip is observed for k 2 > 3 ^1 + jf^yr^-(A V) 2 ^j k 2 
and the equilibrium state arises when equality 
holds. The pattern of amplitude and period of os¬ 
cillation are identical to that for the bubble (Figs. 
5 and 6 ). Figure 5 shows the oscillatory behavior 
of the spike structure for different values of surface 
tension while the dependency of the relative veloc¬ 
ity shear is demonstrated in Fig. 6 . 


and 


2 A 


3(1-A) 



5 1 + A 
16 1-4 


(AV) 2 


[(&) asymp] spike — 0 


provided k 2 < 3 ^1 + pi p ^ pi (AE) 2 ^ k 2 . 


(29) 


IV. Conclusion 

In this paper, we have studied a potential flow 
model to describe the nature of the nonlinear struc¬ 
ture of a two-fluid interface under the combined 
action of Rayleigh Taylor and Kelvin-Helmholtz 
instabilities due to surface tension. The analytic 
expressions for bubble and spike growth rates at 
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Figure 6 : Spike- variation of £i, £2, £3, £4 and £5 
with r. Initial value £1 =-0.1, £2 =0, £3 =0.05, 
£4 = 0, and £5 = 0 with p h = 3, pi = 2, ££=20, 
v h =0, V = 0 (line), 14=0.5, Vj = 0.1 (dot), V h =l, 
Vi = 0.1 (dash), 14=1.5, V[ = 0.1 (dash-dot). 


asymptotic stage are obtained for arbitrary Atwood 
number and velocity shear. Surface tension be¬ 
comes a stabilizing factor of the instability, pro¬ 
vided it is larger than a critical value. In this case, 
oscillatory behavior of motion described by numer¬ 
ical integration of governing equations. The nature 
of oscillations depends on both surface tension and 
relative velocity shear of two fluids. On the other 
hand, below the critical value, surface tension dom¬ 
inates the growth and growth rate of the instability. 
This result is expected to improve the understand¬ 
ing of the stabilization factor for the astrophysical 
instability. 
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